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Prompt heavy quarkonium production in association with a massive 

(anti) bottom quark at the LHC 

Li Gang, Wang ShuangTe, Song Mao, Lin JiPing 
School of Physics and Material Science, Anhui University, Hefei, Anhui 230039, P.R.China 

In this work, we investigate the associated production of prompt heavy quarkonium and a massive 
(anti)bottom quark to leading order in the nonrelativistic factorization formalism at the LHC. We 
present numerical results for the processes involving J/f/j, ^cj, T and Xbj- From our work, we find 
that the production rates of these processes are quite large, and these processes have the potential 
to be detected at the LHC. When pr is smaller than about 10 GeV, the cc^Sq ] state give the main 
contribution to the pr distribution of prompt J/ip with a massive (anti)bottom quark production. 
For the process pp^T + b(b), the contribution of the color-singlet mechanism is larger than that 
in the color-octet mechanism at low pr region. We also investigate the processes pp — > \ c j + b(b) 
and pp — > XbJ + b(b), the pt distributions are dominated by the color-octet Fock state contribution 
at large pr- These processes provide an interesting signature that could be studied at the LHC, 
and the measurement of these processes is useful to test the color-singlet mechanism and color-octet 
mechanism. 

PACS numbers: 12.38.Bx, 13.60.Le, 13.85. Ni, 14.40.Pq 



I. INTRODUCTION 

A heavy quarkonium is made of a heavy quark and a heavy antiquark; its decay and production can be 
factorized into a short-distance part, which can be calculated in QCD perturbatively, and a long-distance 
part, which is governed by nonperturbative QCD dynamics. So the study of heavy quarkonium offers 
a good testing ground for investigating of strong interaction in both perturbative and nonperturbative 
level in high-energy physics. 

In the early days, based on the factorization formalism, color-singlet mechanism(CSM) [1| was used to 
describe the production and decay of heavy quarkonium, and obtained some phenomenological successes. 
However, the CSM has encountered many difficulties in various theoretical [2J and experimental aspects 
[3], such as the appearance of a logarithmic infrared divergence in the case of (next-to-leading-order)NLO 
P-wave decays into light hadrons and the huge discrepancy of the high-pT J/ip production between the 
theoretical prediction and the experimental measurement at the Tevatron. To solve these formal and 
phenomenological problems of CSM, a rigorous theory, nonrelativistic QCD(NRQCD) 4] was proposed 
by Bodwin, Braaten and Lepage(BBL). In the NRQCD, the idea of perturbative factorization is retained, 
the processes of production and decay of heavy quarkonium are also separated into two parts: a short- 
distance part and long-distance matrix elements (LDMEs), which can be extracted from experiments. 
The relative importance of the LDMEs can be estimated by means of the velocity-scaling rules [g| . 

The vital difference between NRQCD and the traditional color-singlet model is that, in NRQCD, the 
complete structure of the quarkonium Fock space has been explicitly considered, and NRQCD predicts the 
existence of color-octet (CO) processes in nature. This CO theory allows that QQ pairs can be produced at 
short distance in CO states and subsequently evolved into physical, color-singlet (CS) quarkonium through 
emission of soft gluons at nonperturbative process @. Being introduced to the color-octet mechanism 
(COM), NRQCD has absorbed the infrared divergences in P-wave 0, 0, @] and D-wave [9|, [l0| decay 
widths of heavy quarkonium and successfully reconciled the orders of magnitude discrepancy between the 
experimental data of J/i/j production at the Tevatron (ll) and the leading order (LO) CSM theoretical 
predictions at large pr- 

Recently, in order to clarify the validity and limitation of the NRQCD formalism, substantial progress 
has been achieved in the calculations of heavy quarkonium production. The DELPHI data favor the 
NRQCD COM prediction for the 77 — > J/tp + X process [12J. Similarly, the recent experimental data 
on the JM> photoproduction of HI [13| are fairly well described by the complete NLO NRQCD cor- 
rections [lj], which give strong support to the existence of the COM. However, at B-factories, a series 
of processes were calculated up to the QCD NLO corrections in the CSM [15l - [l7| . Together with the 
relativistic correction [18|, [l9| , it seems that most experimental data could be understood. Additionally, 
the J /ip polarization in hadroproduction at the Tevatron [2(| and photoproduction at the HERA [21] also 
conflict with the NRQCD predictions. So the existence of the COM is still under question and far from 
being proven. Therefore, the further test for the CSM and COM under NRQCD in heavy quarkonium 
production is still needed. 

In order to investigate the effects of the CSM and COM in heavy quarkonium physics, more processes 
of heavy quarkonium production and decay should be studied. At the LHC, with more attention for the 
heavy quarkonium production, much work has been done in these regions [22j, |2j| , and some channels 
have been calculated to the NLO [24T - [26| . As quarkonium can be identified by using their purely leptonic 
decays [23, [23], the bottom quark can be identified by reconstructing secondary vertices, and the high-py 
bottom quark can be tagged with reasonably high efficiency at the LHC, meanwhile the observation of 
a bottom quark with high-p^ can reduce the backgrounds of the heavy quarkonium production. Thus, 
it is also very interesting to study the prompt heavy quarkonium production associated with a massive 
(anti)bottom quark at the LHC. These processes have the potential to be detected and can provide an 
interesting signature at the LHC. The measurement of these processes is useful to test the CSM and 
COM. 

In this work, we perform the calculations for prompt heavy quarkonium production in association with 
a massive (anti)bottom quark at the LHC. Within this work, the mass of bottom quark is retained in 
all the partonic processes. The paper is organized as follows: We present the details of the calculation 
strategies in Sec. II. The numerical results are given in Sec. III. Finally, a short summary and discussions 
are given. 



II. THE DETAILS OF THE CALCULATION 



The purpose of this paper is to study the associated production of prompt heavy quarkonium and 
a massive (anti)bottom quark to LO in the NRQCD factorization formalism at the LHC. We denote 
the heavy quarkonium as Q. As we know, the cross section for the g(pi) + b(p 2 ) — > Q{k 3 ) + b{ki) 
partonic process in the SM should be the same as that for its charge conjugate subprocess g(pi) + b{p2) -¥ 
2(^3) -\-b(ki) , and the luminosity of the bottom quark in a proton is same as that of the antibottom quark. 
Therefore, the production rates of the Qb and the Qb are identical at the LHC. In the following sections, 
we present only the analytical calculations of the related partonic process g(pi) + b(jpi) — > Q(k s ) + 6(fc4) 
and the parent process pp — > Q + b unless otherwise indicated. 

The cross section for the pp — > Q + b process is expressed as 



a(pp^Q + b) = fdx 1 dx 2 J2(0 Q [n})a(gb^QQ[n}+b) 

J n 

x [G g/A (x 1 ,^ f )G b/B (x 2 ,Hf) + (A <r> B)] 



(1) 



Here a [gb — > QQ[n] + b) describes the short-distance production of a heavy QQ pair in the color, spin and 
angular momentum state n. When Q = c, Q represents charmonium and when Q = b, Q is bottomonium. 
(0^[n\) is the long-distance matrix element, which describes the hadronization of the heavy QQ pair 
into the observable quarkonium state Q. Gi, t g/A,B are the parton distribution functions. A and B refer 
to protons at the LHC. The indices g, b represents the gluon and bottom quark, respectively. 

The short-distance cross section for the production of a QQ pair in a Fock state n, a\gb — > QQ[n] + b], 
is calculated from the amplitudes which are obtained by applying certain projectors onto the usual QCD 
amplitudes for open QQ production. In the notations of Ref. [8| : 
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where A denotes the QCD amplitude with amputated heavy-quark spinors, the lower index q represents 
the momentum of the heavy-quark in the QQ rest frame. Hq/i are spin projectors onto the spin singlet 
and spin triplet states. Ci/ S are color projectors onto the color-singlet and color-octet states. £ a and £ a p 
represent the polarization vector and tensor of the QQ states, respectively. 

Then the LO short-distance cross section for the partonic process g(pi) + b(p2) — > QQ[n\{k^) + ^(fczi) 
is obtained by using the following formula: 



a(gb^QQ[n] + b) = 



1 



16nPN col N v 



d i/ JAjQM 



(3) 



The summation is taken over the spins and colors of initial and final states, and the bar over the summation 
denotes averaging over the spins and colors of initial partons. The Mandelstam variables are defined as: 
s = (pi +P2) 2 , t = (pi — fc.3) 2 , u = (pi —k^) 2 . N co i and N po i refer to the numbers of colors and polarization 
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form in the Appendix. 
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In this paper, we present all contributing partonic cross sections in analytic 






(4) (5) 

FIG. 1: The LO Feynman diagrams for prompt charmonium production in association with a massive bottom 
quark at the LHC. 




(9) (10) 

FIG. 2: The LO Feynman diagrams for prompt bottomonium production in association with a massive bottom 
quark at the LHC. 
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TABLE I: Values of k in the velocity-scaling rule (O [n]) oc v k for the leading QQ Fock states n pertinent to Q. 



In the case of prompt charmonium production in association with a massive bottom quark, there are 
five Feynman diagrams that contribute to this process at LO; we present them in FigQ] There are 



10 Feynman diagrams for prompt bottomonium partonic process, which are drawn in Figj2] In our 
calculations, the QQ Fock states contributing at LO in v for Q are shown in Table [J] 
Following the heavy-quark spin symmetry, these multiplicity relations of LDMEs 

< £>V>("S)p p j8)j > = ( 2 J+ 1) < 0^ nS) [ 3 p( 8) ] >, 

< O r ^[ 3 P^} > = (2J + 1) < T ( nS )[8pW] >, 



< O xc > J fPP] > = (2J+ 1) < O x «°[ 3 P (1) ] >, 



< 



0» J [ 3 S| 8) ] > = (2J+ 1) < 0x«°[ 3 sf >] > (4) 



be assumed satisfied (4J. 



III. NUMERICAL RESULTS 



In our numerical calculations, we focus on the cases Q — J/tp,XcJj Y,Xbj- These quarkoniums can 
be efficiently identified experimentally, and their LDMEs are relatively well constrained |29H33f . The 
?y c , rjb, h c and hb mesons are more difficult to detect experimentally, and we will give the differential cross 
sections relevant to these quarkoniums analytically in our Appendix. 

We take CTEQ6L1 parton distribution functions t 34] with a one-loop running a s in the LO calculations, 
and the corresponding fitted value a s (Mz) = 0.130 is used for our calculations. The factorization scale 

is chosen as \x$ = tot, where tot = \ (PjZ ) + ttvq is the Q transverse mass. The masses of heavy quark 

are set as m c = 1.5 GeV and mb — 4.75 GeV. 

For CO LDME of J/ip, it can be extracted from experimental data or from lattice QCD calculations. 
Before lattice QCD giving out the results, the color-octet matrix elements are determined only by fitting 
the theoretical prediction to experimental data. As done in the literature [13, [30] , the color-octet matrix 
elements are extracted from experimental data of J/tp production in e + e~, ep collisions and hadron 
collisions. In Ref. 19], through the analysis of the process e + e _ — > J/ip + X non - cS at B factories, they 
give a very stringent constraint on the CO contribution, and imply that the values of color-octet matrix 
elements are very small, but their results may not be consistent with the naive velocity-scaling rules. 
Recently, M. Butenschoen and B.A. Kniehl performed a multiprocess fit of the CO LDMEs [30], give 
more global fit results, and their results are consistent with NRQCD scaling rules. In our work, though 
we performed a LO calculation, we will tentatively choose the color-octet matrix elements obtained by 
the M. Butenschoen and B.A. Kniehl's NLO fit results [30] 

<O j ^[ 1 sP]> = (4.76 ± 0.71) x W- 2 GeV 3 , 

<0 J/lp [ 3 s[ 8) }> = (0.265 ± 0.91) x 10- 2 GeV 3 , 

< O j ^[ 3 pP] > = (-1.32 ±0.35) x W- 2 GeV 5 , (5) 

as our input parameters. The relation of the CS matrix elements < O Xc0 [ 3 Pq ] > of Xco with the P- 
wave function at the origin is the formula < O Xc0 [ 3 P (1) ] >= ^|i?p(0)| 2 , and we choose \R' P (0)\ 2 = 
0.075GW 5 from the potential model calculations [3l[. For the CO matrix element < O Xc0 [ 3 S[ S '] >, we 
use < O x "°[ 3 s{ 8) ] >« 2.2 x 10" 3 GeV r3 as our input parameter [Hj. For the NRQCD matrix elements 
of bottomonium, the CS matrix elements are taken from the potential model calculations of [32|, and 
the CO matrix elements are determined from the CDF data (33[. In our calculations, the relation of CS 
matrix elements with the conventions matrix elements of Bodwin-Braaten-Lepage is always considered 
carefully [8J. 

In our numerical calculations, we summed the contributions of pp — > Q + b and pp — > Q+b. If Q = XQJ-> 
we also summed the contributions of XQOiXQi an( i XQ2- We take the constraints of \y \ < 3 for heavy 
quarkoniums, and PTMb) > 5 GeV, | $/&(&) I < 3 for b(b) quark. The colliding energy used in this paper is 
UTeV. 
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FIG. 3: The LO distributions of p T and y '* for the pp —>■ J/ip + 6(6) process, and the contributions of the 
cc^S^] , cc[ 3 S[ 8) ] and cc[ 3 PJ 8) ] Fock states at the LHC. 



In Fig|3l we present the LO distributions of p T and y J ^ for the process pp — > J/tp + b(b) at the LHC. 
At LO, the associated production of prompt J/ip with a (anti)bottom quark is forbidden in the CSM, 
and there is only CO contribution in this process up to the a^v 7 order within the NRQCD framework. 

For comparison, we also depicted the contributions of the cc^Sg ], ccpS^ ] and cc[ 3 Pj ] Fock states in 

these figures. For the LDME of cc[ 3 Pj ] (J = 0, 1, 2) Fock states are negative, the contribution of these 

Fock states are negative too. In FigJ3J we present the px distributions ~da/dp T '^[cc[ 3 Pj ]] in the pt 

distribution figure. From the figure we can see that when pt is smaller than about 10 GeV, the 1 S state 
gives the main contribution to the pr distribution of prompt J/ip with a (anti)bottom quark production. 

With the pt of J '/ip increase, the contributions of cc^Sq ] and cc[ 3 Pj ] Fock states quickly decrease. 

The differential cross section is dominated by the cc[ 3 Sf ] Fock-state contribution at large pr region. In 

the range of 5 GeV < p 1 ^ < 50 GeV, the da/dp 1 ^ is in the range of [0.486, 3276.396]p6/GW, and it 

reaches the maximum when p T = 5 GeV. 



10000 






I 


- 


"1000 






ipl'l 


~: 








tot 


: 


1_ 100 




10 

1 






' 


- 



6 6000 



40 50 

P^(GeV) 



/^ 




1 1 


^. 


,,-''' 




~---.^ 


~ -. 


3 sf 

V™ 






tot 















FIG. 4: The LO distributions of p£ cJ and y XcJ for the pp 
cc[ 3 S[ S) ] and ccfP^] Fock states at the LHC. 



XcJ + b(b) process, and the contributions of the 



The curves for the LO distributions of p^r 7 and y XaJ for the process pp — > XcJ + b(b) at the LHC are 
drawn in Fig|4l From the figure we can see that the contribution from CS is about the same order of 
magnitude with CO at pt — 5 GeV. The CO contribution dominates over production at the large pr 
region, and it decreases much more slowly than that of CS as pr increases. From Fig[3] and FigQ] we 



can find that the cross section of J/tp + 6(6) and \c,j + 6(6) associated production mainly comes from the 
ccpiSj ] Fock state contribution at large px region. The cross section ratios of J/ip + 6(6) and Xcj + 6(6) 
associated production can be estimated approximately by < O j ^[ 3 S[ '] > /J2j=o < XcJ [ 3 S{ '] > 
w 0.134 at large pt region [29|, [30] - So when we consider the process prompt J/tp with a (anti)bottom 
quark associated production at the LHC, the indirect prompt production for this process via radiative 
decays of Xcj should be included. 
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FIG. 5: The LO distributions of p T and y for the process pp 
bbfS^], bb[ 3 S[ 8) ] and bb[ 3 P ( j 8) ] Fock states at the LHC. 
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The LO distributions of pf and y for the process pp — >• T + 6(6), and the contributions of the 



: (8)i 
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,(8) 



bb[ 3 S\% bb^S^ 1 ], 66[ 3 S"} SJ ] and 66[ 3 P) SJ ] Fock states at the LHC are illustrated in FigEJ We can see 



? (i) 



(8)i 



that the differential cross section of bb[ 3 S\ ] and 66[ 3 Sf ] Fock states give the main contribution to the 
distributions of p^ for the process pp — > T + 6(6) at the LO. The 66[ 3 Sf ] state give the main contribution 
to the pt distribution at low px region. With the increase of pj,, the contribution of bb[ 3 S{ ] Fock state 
quickly decrease and the contribution of 66[ 3 S*} ] Fock state becomes more important. The contribution 



J3c(8) 



3c(l)l 



of 66 [ 3 S\ ] Fock states starts to be comparable to 66 [ 3 S\ ] contribution for px — 22 GeV. The differential 
cross section is dominated by the 66[ 3 S'J ] Fock state contribution at large px region. 

In this paper, we have only considered the LO process pp — >• T + 6(6) + X at a 3 . In Ref. [23], the 
authors have evaluated the color-singlet contribution to the process pp — >• T + 66. One of the properties of 
pp — > T + 66 is that its CS contribution (which is at cQ) already includes p^ topologies which dominate 
at large pr- If we consider the process pp — > T + 6(6) + X at af, the process pp — > T + 66 may give 
important CS contribution in large pt region. 

Finally, we give the LO distributions of p^ bJ and y XbJ for the process pp — > Xbj + 6(6), and the 

contributions of the 66[ 3 «S'} ] and 66[ 3 Pj '] Fock states at the LHC in FigJBl From these figures we can 
see that the contribution from CS is larger than CO when pt < 6 GeV. The CO contribution dominates 
over production at the large pt region, and it decreases much slower than that of CS as pt increases. In 
the range of 5 GeV < pf J < 50 GeV, the da/dpf'is in the range of [2.095, 3022.305]p6/Gey, and it 

reaches the maximum when p T = 5 GeV. 



IV. DISCUSSION AND SUMMARY 



In this paper, we investigate the associated production of prompt heavy quarkonium with a massive 
(anti)bottom quark to LO in the NRQCD factorization formalism at the LHC. We have considered all 
experimentally established heavy quarkoniums, with 2S+1 Lj — x 5o, 3 Si, x Pi, 3 Pj, and listed the differential 
cross sections relevant to these quarkonium analytically in our Appendix. We present the numerical 
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FIG. 6: The LO distributions of p^ bJ and y XbJ for the pp — ¥ \bJ + b(b) process, and the contributions of the 
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bb[ d S\ '} and bb[ 6 P) '] Fock states at the LHC. 



predictions of the differential cross sections of the pj, and rapidity y^ for Q — J/ip,XcJ,~£>XbJ at the 
LHC. Wc find that the associated production of prompt Q = J/ip, XcJ, Y, Xbj and a massive (anti)bottom 
quark at the LHC have the potential to be detected. When pt is smaller than about 10 GeV, the cc^Sq ] 
state give the main contribution to the px distribution of prompt J ftp with a massive (anti)bottom quark 
production. For the process pp — >• T + b(b), the contribution of the CSM is larger than that in the 
COM at low pt region. Wc also investigate the processes of pp — > XcJ + b(b) and pp — > Xbj + b(b), 
the Pt distribution are dominated by the CO Fock state contribution at the large pt region. These 
processes provide an interesting signature that could be studied at the LHC, and the measurement of 
these processes is useful to test the CSM and COM. 
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Appendix 

In this Appendix, we list the differential cross sections da/dt for processes gb — > QQ[n] + b. Our results 
read 
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1994s 3 * 5 + 3007s 2 * 6 + \88sV + 8£ 8 )], 



tlgb^bb 
dt l 



lc(8) 
°0 



32Am b (m 2 - s)H 2 (-Aml + i) 2 (-5m 2 + s + t) 4 s 2 



[900000m? 2 - 750m 20 



10 



(5760s + 2569i) + 4m? (2124000s 2 + 1514925s* - 292642?) 
(-8812800s 3 - 6758670s 2 * + 3117666s? + 1530353?) + 8m 



,16 



(649080s 4 + 282870s 3 * - 797714s 2 ? - 401132s? 



72893? 



4m, 



12 



(-440640s 5 + 432465s 4 * + 1876374s 3 ? + 1308895s 2 ? + 231988s? 
+29447?) + 4m?(84960s 6 - 533250s 5 * - 1391312s 4 ? - 1158952s 3 ? 
-439340s 2 ? - 33147s? - 3587?) + s 2 t(s + *) 2 (90s 6 + 270s 5 * 
+363s 4 ? + 276s 3 ? + 125s 2 ? + 32s? + 4?) + m 8 /-34560s 7 
+976500s 6 * + 2600172s 5 ? + 2623254s 4 ? + 1433016s 3 ? + 381832s 2 ? 
4mg(360s 8 - 60660s 7 * - 186316s 6 ? 
s 3 ? - 13166s 2 ? 

+2m 4 *(17325s 8 + 63180s 7 * + 97678s 6 ? + 87168s 5 ? + 50438s 4 ? 
+17006?? + 1993s 2 ? + 32s? + 2?) - 4m 2 s*(675s 8 + 2910s 7 * 
+5448s 6 ? + 5912s 5 ? + 4167s 4 ? + 1915s 3 ? + 484s 2 ? + 39s? + 2?)], 



+9382s? + 1141? 



-233112s 5 ? - 165986s 4 ? - 72110s 3 ? - 13166s 2 ? - 78s? - 19? 



da 
di 



(gb 



bb 



3 5 (D 



-32a 3 7r 2 
243m 6 (m 2 - iW^2 

+u) 
+23389*?) 



'■b - «?( m 
2s 2 u 2 {s 



u) 4 (-2m 2 



+ u) 3 (9s 2 



s + u 

.-.2 



2^2 [5466m? - 18668m?(s 



+m 
-m 8 



■1° (4609s 4 



., , „, ,.,., + 16sw + 9u 2 )+m?(23389s 2 + 62070sM 

m?(14251s 3 + 73585s 2 u + 73585si? + 14251w 3 ) 
:•>.-,, :;.-. + 7458832^2 + 432593^ + ^qqq^) 

35184s 3 t? + 35184s 2 w 3 ' 

2j»/nn«fi , nnn-5- , wr> -1-4-2 , nn A 



*(807s 5 + 14717s 4 



u - 



14717su 
+807i?) + m 2 su(20s 6 + 228s 5 w + 701s 4 w 2 + 984s 3 i? 
+701s 2 u 4 + 228si? + 20i?) + mg(72s 6 + 3179s 5 



u 
|-12444s 3 u 3 

M753s 5 i? 



v, + 9455s 4 u 2 
a 4 (2s 7 



406s 6 u 



f 9455s 2 u 4 + 3179su° + 72m") - m% x - 
3043s 4 i? + 3043s 3 m 4 + 1753s 2 w 5 + 406st? + 2u 7 )], 



da 
di 



gb — > bb 



3 5 (8) 



[66174m 20 - 226934m? 



1944m 3 Xm 2 - s) 4 (ml - u) 4 (-2m 2 + s + u) 2 s 2 
(s + u) + m? (260287s 2 + 806952su + 260287i?) - 10m?(10975s 3 
+96181s 2 *i + 96181su 2 + 10975t?) + 9?i?(4s 4 - s 3 u + 8s 2 u 2 - si 
+4w 4 ) + 2m?(6155s 4 + 216310s 3 u + 584812s 2 t? + 216310si? 
+6155u 4 ) - 2m^s 2 i?(90s 5 + 493s 4 w + 468?i? + 468s 2 u 3 + 493sw 



-90u 



r,"n 



2m?(1548s 5 + 27601s 4 u + 266181s 3 t? + 266181s 2 u 3 



-27601su 4 



1548u 5 ) + 2m 4 su(94s 6 + 534s 5 w + 3748s 4 u 2 - 477s 3 w 3 
+3748J 2 ?? + 534su 5 + 94t?) + mg(285s 6 + 9920s 5 u + 60701s 4 w 2 
+245736?i? + 60701s 2 w 4 + 9920si? + 285i?) - 2m^(22s 7 + 563s 6 u 
+4424s 5 t? + 13715s 4 ?? + 13715s 3 u 4 + 4424s 2 i? + 563su 6 + 22u 7 )], 



dt \ 



bb 



^ 



(i) 



b = 



64c?t? 



243(m 2 -s) 5 (-5m 2 + s 
(10071936s + 1099435*) 
(-6754944s 2 + 28431108s* + 367984? 



[2227744r 



-2m( 



-*) 5 (4m 3 -m b t) 3 s 2 
s 3 t 4 (s + *) 3 (2s 2 + 2s* + 9?) + m 20 



,18 



(12743424s 3 



11 



-25261690s 2 * - 16277654s? + 488169?) - 2m£ 6 (3126432* 



+13119432s J * + 14986264sV - 2170280s* J + 183541* 



4to 



14 



(409536s 5 + 2567733s 4 * + 5640894s 3 ? + 4564301s 2 ? - 49258s? 
+32054?) + m 2 b s 2 t(s + *) 2 (2s 6 + 6s 5 * + 35s 4 ? + 108s 3 ? + 249s 2 ? 
+334s? + 27i 6 ) - to?(257664s 6 + 2250792s 5 * + 7210544s 4 ? 
+10902208s 3 ? + 6680880s 2 ? + 222148s? + 27053?) + 2to? 
(11904s 7 + 152098s 6 * + 651630s 5 ? + 1443207s 4 ? + 1656100s 3 ? 
+780410s 2 ? + 36206s? + 1776?) + mg£(1474s 8 + 10168s 7 * 
+36644s 6 ? + 87208s 5 ? + 125184s 4 ? + 88392s 3 ? + 22168s 2 ? 
+836s? + 9?) - to 4 s*(60s 8 + 448s 7 * + 1888s 6 ? + 5568s 5 ? 
+11620s 4 ? + 14011s 3 * 5 + 7614s 2 * 6 + 1182s* 7 + 27?) - m s b (992s 8 
+26384s 7 * + 143952s 6 ? + 424880s 5 ? + 736548s 4 ? + 663440s 3 ? 
+235347s 2 ? + 10823s? + 269?)], 



da 
di 



%(***['* 



>(8) 



b) = 



[32400000m 28 - 5400m 26 



972(to 2 - s) 5 ?(-5to 2 + s + *) 5 (4mf - m b t) 3 s 2 

(36000s + 57209*) + 2m 24 (249480000s 2 + 561226320s* + 23814083?) 

- 2mf (357696000s 3 + 843783912s 2 * + 267419628s? - 77026763?) 
+ m 20 (628819200s 4 + 1367173728s 3 * + 912846660s 2 ? - 48175560s? 

- 99056939?) 4 

+ 60998614s 2 ? 

+ 69346368s 5 * - 

- 37003547s? - 

S3+3 



m? (-351226368s 5 - 604100808s 4 * 



89194312s? + 28705941?) 



m 



i(i 



504111528s 3 ? 

is6 



(125763840s 6 



243821166s 4 ? 



449357484s 3 ? 

r.3i2/ s , ,T\3/ifioa6 , aqo^} , ,.„.,;.l,'2 



283946945s 2 ? 



5575172?) - s A t\s + *) 3 (162s 6 + 486s 5 * + 693s 4 ? 



576s 3 ? + 323s 2 ? + 116s? + 36?) + 4to?(-7153920s 7 + 20626056s 6 * 
112036932s 5 ? + 207447531s 4 ? + 131086836s 3 ? + 43878881s 2 ? 
3359759s? + 243128?) + to? (3991680s 8 - 55432512s 7 * 
267148776s 6 ? - 604619712s 5 ? - 587249390s 4 ? - 245568940s 3 ? 
61736796s 2 ? - 3743122s? - 151151?) + m 2 b s 2 t(s + *) 2 (648s 8 
8424s 7 * + 29366s 6 ? + 50244s 5 ? + 47759s 4 ? + 26802s 3 ? - 



s 2 ? 



,10/ 



91385904s 7 ? 



+ 1840s*' + 108?) + TO?(-311040s y + 17461656s * 

+ 243458012s 6 ? + 322380960s 5 ? + 207346758s 4 ? + 66243932s 3 ? 

+ 13294060s 2 ? + 642412s? + 16665?) + to 8 , (10368s 10 - 3252960s 9 * 

- 19662902s 8 ? - 60570920s 7 ? - 99701250s 6 ? - 88396106s 5 ? 

- 42116100s 4 ? - 10893998s 3 ? - 1728483s 2 ? - 64271s? - 1112?°) 
+ 2TO 6 i(183384s 10 + 1338660s 9 * + 4750079s 8 ? + 9250520s 7 ? 

+ 10317034s 6 ? + 6712106s 5 ? + 2525028s 4 ? + 540096s 3 ? + 65441s 2 ? 



IS8X2 



1852s? + 18? u ) - to 4 s*(23328s 1u + 220572s 9 * + 904920s 8 * 



2024407s 7 ? 
62406s 2 ? 



2686780s 6 ? 



2197140s 5 ? 



1116238s 4 ? 



346213s 3 ? 



5448s? + 108?°)], 



^(gb- 
dt l 



bb 



3p(l) 



b) = 



4a 3 7? 



81 (to 2 - s) 5 ?(-5m 2 + s + *) 5 (4to| - to?) 3 s 2 



[25200000to 28 - 600to 26 



12 



(252000s + 346111*) 

+ 18ml 4 (21560000s 2 



3s 3 i 2 (s + *) 3 (3s 2 + 3s*- 
f 26837680s* + 6879209? 
58438984s 2 * + 132101836s* 2 - 31650279* 3 ) + 
52223904s 3 * - 596240052s 2 * 2 



4* 2 ) 2 (2s 2 + 2s* + * 2 



,20 



6m 22 (92736000s 
(489081600s 4 



::'■'> 



3ml 8 (-91058688s 5 + 34245912s 4 * 
29438668s* 4 + 71917105* 5 ) 



475090104s* 3 - 316136919* 4 ) 

543856590s 2 * 3 



252722648s 3 * 2 



,16 



(97816320s 6 - 160017984s 5 * 



491228090s 4 ? - 877626916s 3 * 3 - 1687844825s 2 * 4 - 268763655s* 5 
87840460* 6 ) + 4to^ 4 (-5564160s 7 + 30448824s 6 * + 75718764s 5 * 2 
140914224s 3 * 4 + 246581099s 2 * 5 + 40404191s* 6 



66864981s 4 * 3 



151809336s 6 * 2 



5840348*0 + m? (3104640s* - 54127296s'* 
114986592s 5 * 3 - 65176694s 4 * 4 - 238851724s 3 * 5 - 359985180s 2 * 6 
52213906s* 7 - 4131047* 8 ) + m 2 b s 2 i{s + *) 2 (504s 8 + 3960s 7 * 
8522s 6 * 2 + 3336s 5 * 3 - 4911s 4 * 4 - 1942s 3 * 5 + 791s 2 * 6 + 416s* 7 



A 7+2 



1552* s ) + mn-241920s J + 14973288s 8 * + 51568272s'* 2 
52527524s 6 * 3 + 13171944s 5 * 4 + 8800346s 4 * 5 + 67284444s 3 * 6 
85100380s 2 * 7 + 10250336s* 8 + 471419* 9 ) + to 8 , (8064s 10 
2633760s 9 *- 11282354s 8 *" 2 - 15962360s 7 * 3 - 7053666s 6 * 4 
1332094s 5 * 5 - 874708s 4 * 6 - 12363406s 3 * 7 - 13033959s 2 * 8 
1228903s* 9 - 31560* 10 ) + 2mgt(144360s 10 + 773580s 9 * 
1467393s 8 * 2 + 1038192s 7 * 3 + 56862s 6 * 4 - 231222s 5 * 5 + 62812s 4 * 6 
716948s 3 * 7 + 623079s 2 * 8 + 41600s* 9 + 472* 10 ) - m^s*(18144s 10 
124628s 9 * + 311800s 8 * 2 + 318915s 7 * 3 + 73212s 6 * 4 - 68716s 5 * 5 
37362s 4 * 6 + 24645s 3 * 7 + 98510s 2 * 8 + 67280s* 9 + 2448* 10 )], 



da 
di 



gb — > bb 



3p(8) 



b) = 



a 3 7r 2 



324(m 2 



5to| 



s)H 2 { 
21000to 26 (36000s 



[126000000m 28 



+ s + i) 5 (4m 3 — TO;,*) 3 s 2 

5111*) + 6m 24 (323400000s 2 - 10772400s* 



179385097* 2 ) - 6m 22 (463680000s 
199413487* 3 ) 4 
765731304s* 3 - 
839359544s 3 * 2 



S3 



66541160s 2 * - 356589988s* 2 
3m 20 (815136000s 4 + 161689440s 3 * - 713740812s 2 * 2 
210248369* 4 ) - 3to^ 8 (455293440s 5 - 60499320s 4 * 
- 998696926s 2 * 3 - 412742996s* 4 - 69091505* 5 



m 



16 



(489081600s 6 - 695580480s 5 * - 3191554274s 4 * 2 - 3864588148s 3 * 3 
- 2108766941s 2 * 4 - 451377411s* 5 - 49514524* 6 ) - 3s 3 * 2 (s + *) 3 (9 
+ 270s 5 * + 363s 4 * 2 + 276s 3 * 3 + 125s 2 * 4 + 32s* 5 + 4* 6 ) + 4to^ 4 
(-27820800s 7 + 147889560s 6 * + 590044860s 5 * 2 + 991176813s 4 * 
+ 673013160s 3 * 4 + 228128459s 2 * 5 + 31137503s* 6 + 2314418* 7 ) 



it:; 



; (15523200s 8 



268977600s 7 * - 
2341866062s 4 * 4 - 1031915260s 3 * 5 - 
1320599* 8 ) + m 2 b sH(s + *) 2 (2520s 8 
68064s 5 * 3 + 64785s 4 * 4 + 48194s 3 * 5 
to^°(-1209600s 
1006700472s 5 * 4 



1077142104s 6 * 2 - 2165379264s 5 * 3 
- 254801556s 2 * 6 - 25728466s* 7 
+ 19800s 7 * + 51098s 6 * 2 
+ 17783s 2 * 6 + 344s* 7 + 388* 8 ) 
'' '- 74797320s 8 * + 324287952s 7 * 2 + 702527876s 6 * 3 
755209754s 4 * 5 + 232438668s 3 * 6 + 46247836s 2 * 7 
3738896s* 8 + 131591* 9 ) + m% (40320s 10 - 13168800s 9 * - 64937498s 8 * 2 



13 



- 149484440s 7 * 3 - 240269802s 6 * 4 - 255755594s 5 * 5 - 139606876s 4 * 6 

- 30968470s 3 * 7 - 5385543s 2 * 8 - 359983s* 9 - 8100* 10 ) + 2m 6 ,* 
(721800s 10 + 4187100s 9 * + 10660281s 8 * 2 + 17776032s 7 ? + 22160286s 6 * 4 
+ 17930202s 5 * 5 + 7326376s 4 * 6 + 1157204s 3 * 7 + 196287s 2 * 8 + 10640s* 9 
+ 118* 10 ) - m 4 s*(90720s 10 + 644420s 9 * + 1920040s 8 * 2 + 3393519s 7 * 3 

+ 4354860s 6 * 4 + 4215260s 5 * 5 + 2625438s 4 * 6 + 805269s 3 * 7 + 84038s 2 * 8 
+ 17240s* 9 + 612* 10 )1. 
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